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ABSTRACT 

Stellar population parameters (e.g. age, metallicity, and stellar abundance ratios) derived from 
spectral line-strengths provide a powerful probe of galaxy properties and formation histo- 
ries. We implement the machinery for extracting line-strengths and 'single-stellar-population- 
equivalent' (SSP-equivalent) stellar population parameters from synthetic spectra generated 
by a hierarchical galaxy formation model. Our goals are (1) to test the consistency of these 
line-strength-derived stellar population parameters with more physically relevant light- and 
mass-weighted parameters for complex, cosmologically-motivated star-formation histories, 
(2) to interpret line-strength observations for early-type galaxies within the context of hierar- 
chical structure formation, and (3) to test the galaxy formation models using stellar popula- 
tion parameters derived from some of the best available samples of observed line-strengths. 
We find that the SSP-equivalent age is related to the light-weighted age in a complicated 
fashion that reflects the influence of recently-formed stars and is poorly correlated with the 
mass-weighted age. We find that the tendency for SSP-equivalent ages to be biased young 
means that 'archaeological downsizing' overstates the 'true', mass-weighted downsizing in 
age with mass. We find however that the SSP-equivalent metallicity closely tracks the mass- 
and light-weighted metallicities, so that observed mass-metallicity relations for old galaxies 
closely reflect the underlying trends. We construct mock catalogues of early-type galaxies in 
a Coma cluster-sized halo and compare them directly to observations of early-type galaxies 
in the Coma cluster The similarity of the SSP-equivalent ages in the observational samples 
and the mock catalogues gives us confidence that the star-formation quenching implemented 
in the hierarchical galaxy formation model produces roughly the correct amount of recent 
star formation. Unfortunately, the current observational samples are either too small or have 
too low signal-to-noise to accurately determine detailed star-formation histories. However, 
the data show that the model has deficiencies: the SSP-equivalent metallicities are too low 
and have the wrong slope as a function of velocity dispersion, and the SSP-equivalent ages 
of the model galaxies may have an incorrect slope as a function of velocity dispersion. These 
problems are indicative both of the simplified chemical evolution prescription currently im- 
plemented in the galaxy formation model and that the star-formation histories resulting from 
the model are incorrect in detail. 

Key words: galaxies: formation — galaxies: evolution — galaxies: stellar content — galax- 
ies: ellipticals and lenticulars — galaxies: clusters: individual (Coma) 



1 INTRODUCTION 

The optical colours of old stellar populations are determined by an 
interplay between the colours of the red giant branch and the main- 
sequence turn-off, which together produce nearly all of the light at 
these wavelengths. This interplay results in the well-k nown age- 
metallicity degeneracy in old stellar populations (e.g.. lO'Cormelll 
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Il986l) : a change by a factor of two in age looks the same in optical 
colours (an d metal-line str engths) as a factor of three change in 
metallicity ( IWorthevl[T994l) . This age-metallicity degeneracy long 
limited - and to some extent, still limits - our understanding of the 
stellar populations of 'red, dead' galaxies like early-type galaxies. 

In a pioneering paper, iRabid ( Il982h demonstrated that it is 
possible to break the age-metallicity degeneracy by using the 
strength of the Balmer lines of hydrogen as a function of a metal 
line. Ever since this breakthrough, the potential power of 'fossil' 
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star formation evidence in local galaxies has been expected to yield 
insights on the formation of these objects comple mentary to direc t 
look-back studies (see, e.g., the recent review bv lRenzin ?'2006^. 
If we could measure the mass-weighted ages of local galaxies, 
we would know when the bulk of their stars were formed, a key 
constraint on galaxy formation models. This approach has yielded 
very interesting results when authors have applied modem stellar 
population models (by, e.g., |B uzzoni et al. 1994; Worthev 1994'; 
Vazdekis et alJl996l ;lBruzual & Chariot 2003; Thomas et al. 2003 



; et all 

SchiavonlbOOTl and others) to high-quality spectral observations of 



early-type galaxies. For example, these studies have demonstrated 
the wide range of galaxy ages (Gonzalez 1993), the dependence 
of stellar population parameters (ag e, metallicity, and abundance 
ratio) on galaxy velocity dispersion jTrager et alj|2000bi) , the de- 
pendence of age and star-formation time-scale on mass (often re- 
ferred to as 'fossil' or 'a rchaeological downsizing': iThomas et al.l 
I2OO5I ; iNelan" et alj|2005l) , and the dependenc e of galaxy age and 



star formation time-scale on environment 



iBemardi et"Zll2006l; Isanchez-Blazguez et al.H2006bl, but 



tradicatory evidence in lThomas et al .II2OO7I) 



Thomas et all l2005l: 



But as age-dating techniques become more routine and con- 
front larger data sets, the pitfalls of this approach become more 
apparent. The most dangerous pitfall is that the easiest and there- 
fore most common application of these techniques assumes that 
galaxies can be meaningfully parametrized by a single 'age'. 
This assumption is in clear violation of the modern picture of 
galaxy formation, in which (early-type) galaxies do not form in 
a single burst, but are built up through a sequence of merg- 



ers (e.g.. 'Toomre 1977: White & Rees 


1978 


: Blumenthal et al.l 


1984: Kauffmann et alj 19931: Cole et al 


1994 


: Kauffmann! 1994 


Kauffmann & CharlotI 1 19981: iDe Lucia et al. 2006. among manv 



others). As shown schematically bylTrager et al. ( 2 000h) and quan 



others). As shown schematically by 
titatively bv lSerra & Tragej feOOTn . 



the addition of a small fraction 
of young stars to an old population strongly biases the apparent age 
of a galaxy. As a simple example, the addition of 2 per cent by mass 
of 1 Gyr-old s tars to a 12 G yr-old population results in an apparent 
age of 5 Gyr cTrager et al.ii2008 ). The extreme sensitivity of stellar 
population ages to recent star formation leads to ambuigities in in- 
terpretation: is a galaxy truly young or just composite, as would be 
expected from our current understanding of galaxy formation'^ 

Furthermore, accurate age dating, even assuming a single stel- 
lar population, requires very accurate line strengths and therefore 
very high signal-to -noise spectra (S/N > 100/A : iTrageil [19971 ; 
ICardielet alii 19981 l2003l : iKuntschner 67^1200 ih . This limits its 
usefu lness so far eith er to small data s e ts (e.g., Gon zalez 1993i; 
Kuntschner 2000; Thon ias et al.l l2005l : ISanchez-Blazguez et al.l 
2006a ; .Trager et al. ,200^ or to the need to combine dozens-to- 



hundreds of galaxies together, thus losing inf ormation about dis- 
persions in stellar p opulation properties (e.g., Bernardi et al.||200^ ; 
[Graves et al. I l2007l) . It is therefore clear that guidance from simu- 
lations could be very useful for understanding the impact of both 
complex star-formation histories and noisy data on the interpreta- 
tion of stellar population ages, if simulation results could be ana- 



^ Although more sophisticated methods of disenta ngling star-formation 
histories from spectra are becoming available (e.g., iHeavens et alj |2000| ; 
iReichardt etai]|200ll ; lOcvirk et al1l2006l ; iToieiro et alj|2007l) . thevappear 
either to require very hig h signal-to-noise and high-resolution s pectra over 
a long spectral baseline jOcvirk et alj|200 j ; IXoieiro et alj|200"7l) or to pro- 
duce results that are useful on averag e only for very large samples (e.g. 
IPanter et alj2003l ; |jimenez et alj20"o% . 



lyzed, to the greatest extent possible, in exactly the same way as the 
observations. 

We turn to the latest generation of hierarchical galaxy for- 
mation models to provide this guidance. This new generation of 
models implements feedback from active galaxy nuclei, most im- 
portantly in the form of heating from radio jets, to solve the over- 
cooling (e.g., Kauffmann et al. 1993) and star-formation quenching 
(e.g., Som erville et al. 2004) proble ms of the previo us generation 
of models jBower et al. 2006 ; Cattaneo et al.1 12OO6I : ICroton et al] 
l2006l : ISomerville et alJl2008bE . These models successfully predict 
many of the properties of local galaxies, in particular their luminos- 
ity fu nctions, the dichotomy of galaxy colours (e.g.. lStrateva et all 
l200lh . the morphology of the colour-magnitude diagram, and even 
its evolution to 2: = 1. However, the authors of these papers make 
predictions almost solely in colour-magnitude space, which is at 
best a blunt tool: as discussed above, the age-metallicity degener- 
acy makes interpreting the colours of old populations problematic. 

Now that the hierarchical galaxy formation models are doing 
a reasonable job at predicting basic local galaxy properties, we in- 
vestigate using these models to produce 'mock catalogues' to guide 
our interpretation of line-strength observations of local galaxies. 
We further want to subject the models themselves to a sharper test 
than the blunt tool of colours alone. We have two desires in this 
paper to make these goals concrete: (1) to use a well-defined sam- 
ple of galaxies drawn from a well-understood environment that can 
be simply modelled in the context of the hierarchical galaxy for- 
mation models and (2) to bring the models as close to the 'obser- 
vational plane' as possible. To meet the first desire, we use line- 
strength observations from recent samples of early-type galaxies 
in the Coma cluster. We choose the Coma cluster because it is a 
rich, well-studied cluste r whose dark matter halo properties are rel - 
atively well-understood fcokas & Mamonll2003l ; [Kubo et alj2007h . 
The Coma cluster is close enough to allow detailed spectroscopic 
study of a large number of early-type galaxies and yet is far enough 
that the galaxies can be treated almost as 'point sources' (but not 
quite: see the Appendix) in a single dark matter halo, a convenience 
when comparing to the hierarchical models. We will find below 
that the current line-strength observation samples of Coma clus- 
ter galaxies do not provide the ideal galaxy sample for our needs, 
but they are the best we have at the moment. To meet the second 
desire, we create synthetic spectra from the hierarchical galaxy for- 
mation models, based on the predicted star formation and chem- 
ical enrichment histories convolved with stellar population mod- 
els, and run them through the same machinery used to extract line- 
strengths, SSP-equivalent parameters, and error-bars on these pa- 
rameters from the observed spectra. We apply this machinery to 
mock catalogues drawn from simulated Coma-sized halos such that 
we reproduce the Velocity Dispersion Function (VDF) of a partic- 
ular observational sample. These predictions give us the sharp tool 
we need to understand the properties of local galaxies and to test 
the models themselves. 

We have three goals in the current paper: (1) understand- 
ing and calibrating how SSP-equivalent parameters extracted from 
spectral line- widths represent the underlying mass- and luminosity- 
weighted ages and metallicities for galaxies with complex, cosmo- 
logically motivated star formation and enrichment histories, (2) in- 
terpreting line-strength observations of early-type galaxies within 
the context of hierarchical models, and (3) using the observed line- 
strengths for the best available samples of Coma cluster galaxies to 
test the physics in the galaxy formation models. The paper is orga- 
nized as follows. In Section|2]we review the ingredients of the hier- 
archical galaxy formation models, their calibration, and the stellar 
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population modelling. In Section[3]we examine the SSP properties 
of composite populations from synthetic spectra created via the hi- 
erarchical galaxy formation models. In Section|4]we compare mod- 
els of galaxies in the Coma cluster with observations and suggest 
improvements for both models and observations. Finally in Sec- 
tion|5]we summarise our findings. An appendix details the aperture 
corrections needed to compare the resolved galaxy observations to 
the unresolved model galaxies. 



2 HIERARCHICAL GALAXY FORMATION MODELS 

We make u se of a se mi-analyti c galaxy formation model de s cribed 
in detail inlSomerville & PrimacM h999l) . ISomerville et alj j200lh 



and lSomerville et al.l i2008bl . hereafter SOS). 



2.1 Merger trees and substructure 

The models are based on Monte Carlo realizations of dark matter 
halo merger histories, constructed using a slightly modified ver- 
sion of the method of ISomerville & Kolatti jl999h . as described in 
SOS. The merging of sub-halos within virialized dark matter halos 
is modelled by computing the time required for the sub-halo to lose 
its angular momentum via dynamical friction and fall to the center 
of the parent halo. Our treatment takes into account the mass loss 
of the orbiting dark matter halos due to tidal stripping, as well as 
tidal destruction of satellite galaxies. 



2.2 Gas cooling and star formation 

When gas cools, it is assumed to initially settle into a thin exponen- 
tial disc, supported by its angular momentum. Given the halo's con- 
centration parameter, spin parameter, and the fraction of baryons in 
the disc, we can use angular momentum conservation arguments 
to compute the scale rad ius of the exponential disc after collapse 
dSomerville etai]|2008 j) . 

The models make use of fairly standard recipes for the sup- 
pression of gas infall due to a photo-ionizing background, cooling, 
star formation, and supernova feedback, as described in SOS. Qui- 
escent star formation, in isolated discs, is modelled using the em- 
pirical Kennicutt star-formation law CKennicutt., 1 9S9, . 1 99Si) : 



SsFR = A 



Kcnn -^gas 



(1) 



with j4Kcnn = 8.33 X 10^"" (a factor of two lower than the canon- 
ical value, see SOS), Nk = 1.4, Egas the surface density of cold 
gas in the disc (in units of Mq pc~^), and where Ssfr has units of 
Mq yr~^ kpc~^. We also adopt a critical surface density threshold 
Ecrit , and assume that only gas lying at surface densities above this 
value is available for star formation. 



2.3 Galaxy mergers and morphologies 

Galaxy mergers may induce a burst of star formation and destroy 
any pre-existing disc component, depending on the mass ratio of 
the smaller to the larger galaxy, /i. We parameterize the strength 
and timescale of these bursts according to the results of a large suite 
of hydrodynamic simu lations of galaxy mergers (Robertson et al. 
l2004ICoxetal]|2008l) as described in SOS. The burst efficiency 
is a function of the mass ratio ^ and the bulge-to-total fraction of 
the primary galaxy (larger mass ratios yield more dramatic bursts, 
while the presence of a large bulge further suppresses bursts in low 



mass-ratio mergers); the burst timescale is primarily a function of 
the primary progenitor galaxy's circular velocity (lower circular ve- 
locity galaxies produce longer bursts). 

Mergers can also heat and thicken, or even destroy, a pre- 
existing disc component, driving galaxies towards morphologically 
earlier Hubble types. We assume that the fraction of the pre-existing 
stars that is transferred from a disc to a spheroidal component is a 
strongly increasing function of the merger mass ratio /i, such that 
minor mergers with /i < 0.2 have little effect, and major mergers 
with fj, > 0.25 leave behind a spheroid-dominated remnant. 

Morphologies of the model galaxies are determined using 
the bulge-to-tot al i3-band ligh t ratios B / T following the scheme 
described in iKauffmann et al.l (|l993) and Somerville & Primac. 
Il999.) , based on the results of ISimien & de Vaucouleurs (19Si 
Galaxies with B /T < 0.405 are considered to be spirals, galaxies 
with 0.405 < B/T < 0.603 are considered to be SOs, and galax- 
ies with with B/T > 0.603 are considered to be ellipticals. In the 
current paper we consider galaxies with type SO or E to be early- 
type galaxies. Note that we conservatively use B magnitudes un- 
corrected for dust extinction in this calculation. This choice moves 
galaxies to somewhat later morphological types (smaller B /T) as 
the discs are extinguished more than the bulges, so the unextin- 
guished discs are brighter than they would be observed. 

It is important to note here that although the semi-analytic 
model provides predictions of size estimates for disc galaxies 
(see SOS), the model currently does not provide size estimates 
for bulges. Therefore when comparing properties of local bulge- 
dominated (early-type) galaxies to model predictions, we must be 
aware of the possible effects of spatial gradients on these results. 
We discuss this in more detail in Sec.|4]and in the Appendix. 

2.4 Black hole growth, supemovae, and feedback 

In our model, mergers also trigger the accretion of gas onto super- 
massive black holes in galactic nuclei. Each top-level halo in our 
merger trees is seeded with a black hole of mass M^ccd — IOOM0 . 
Following a merger, the black hole is allowed to grow until the 
radiative energy being emitted by the AGN becomes sufficient 
to halt further accretion. This self-regulated treatment of black 
hole growth is based on hydrodynamic simulations including BH 
growth and feedbac k dSpringel et Zll2005l : iDi Matteo etal]|2005l : 
iHopkins etai]|2007l) and is described in more detail in SOS. Be- 
cause the potential well depth of the spheroid determines how large 
the black hole has to grow in order to halt further accretion, this 
model successfully reproduces the black hole mass vs. spheroid 
mass relation at z ~ 0. Energy radiated by black holes during 
this 'bright', quasar-like mode can also drive galactic-scale winds, 
clearing cold gas from the post-merger remnants (see SOS). 

In addition to the rapid growth of BH in the merger-fueled, 
radiatively efficient 'bright mode', we assume that BH also experi- 
ence a low-Eddington-ratio, radiatively inefficient mode of growth 
associated with efficient production of radio jets that can heat gas 
in a quasi-hydrostatic hot halo. The accretion rate in this phase is 
modelled assuming Bondi accretion u sing the isothermal cooling 
flow solution of .Nulsen & FabianI ( l2000i) . We then assume that the 
energy that effectively couples to and heats the hot gas is given by 
i'hcat ~ Khcat?7nT-radioC^, whcrc rfiradio IS the accretion rate onto 
the BH, r; = 0.1 is the assumed conversion efficiency of rest-mass 
into energy, and ACheat is a free parameter of order unity. 

Cold gas may be reheated and ejected from the galaxy, and 
possibly from the halo, by supernova feedback. The rate of reheat- 
ing of cold gas is given by 
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= eo I y I "I., (2) 

where Q^h = 2 is assumed for energy- driven winds (SOSfl The 
heated gas is either trapped within the potential well of the dark 
matter halo and deposited in the 'hot gas' reservoir, or is ejected 
from the halo into the 'diffuse' Intergalactic Medium (IGM). The 
fraction of reheated gas that is ejected from the halo is given by: 



/cjcct(Kir) = [1.0 + (Vvir/Kjcct)""^" 



(3) 



where aojoct = 6 and Vcjcct is a free parameter in the range ~ 100 
-150 km s~^. This ejected gas is reincorporated into the halo on a 
timescale proportional to the halo dynamical timescale (see SOS). 

2.5 Chemical evolution 

Each generation of stars produces new heavy elements. The mass 
of metals produced is dMz = ydm^,, where dm, is the mass of 
stars produced and y is the yield. New metals are added directly to 
the cold gas. When the cold gas is reheated or ejected by supernova 
feedback, the metals are transferred to the hot or ejected compo- 
nent, respectively, in the same manner as the gas. Ejected metals 
are also 're-accreted' by the halos as described in SOS. Note that 
the simple chemical enrichment model used here includes only en- 
richment by Type II supernovae. Therefore the metallicities quoted 
should be interpreted as representing a-type elements. In this study 
we assume that [E/Fe] = 0, that is, that the galaxies have solar 
abundance ratios, and that the metallicities quoted are total metal- 
licities. This assumption has some impact on the results discussed 
below, in particular on the inferred metallicities of the model galax- 
ies. In a forthcoming paper (Arrigoni et al., in prep.) we relax this 
assumption and include a full galactic chemical evolution (GCE) 
model, tracing up to 19 individual elements. We point out below 
where this full GCE model differs from the present models. 

2.6 Calibration 

The models contain a number of free parameters, w hich are sum- 
mariz ed in Table 1 of SOS. We adopt the 'WMAP3' jSpergel et al.l 
I2OO7I) values of the cosmological parameters (SOS). We then set the 
free parameters by matching certain observations of nearby galax- 
ies, as described in detail in SOS. For example, the star formation 
parameters AKenn and Ecrit are constrained by the observed gas 
fractions of nearby spiral galaxies. The supernova feedback param- 
eters are set in order to reproduce the stellar mass fractions in low- 
mass halos ( < IO^^Mq) such that the low-mass or low-luminosity 
slope of the observed stellar mass or luminosity function is repro- 
duced. The yield y is set such that spiral galaxies in Milky Way 
sized halos (2 x IO^^Mq) have approximately solar metallicity. 
The parameter controlling radio mode AGN feedback is set to give 
the minimum amount of feedback required to solve the overcool- 
ing problem and bring the high-mass or high-luminosity end of 
the galaxy stellar mass or luminosity function into agreement with 
observations. The resulting models, which are tuned solely using 
z ^ observations, are also in fairly good agreement with direct 
probes of the mass assembly and star-formation history of the uni- 
verse from look-back studies based on deep cosmological surveys 
(SOS). 



^ Note that there is a typographical error in this equation in SOS. The ex- 
pression given here is the correct one. 



2.7 Stellar kinematics 

To compare model galaxy properties with observations, a use- 
ful observational tracer is the one-dimensional line-of-sight ve- 
locity dispersion, a (hereafter referred to simply as the velocity 
dispersion). Velocity dispersion appears to be the controlling pa- 
rameter of stellar population ob servables, including colour and 
lumi nosity teem ardi et al.M2005!). metallicity and abundance ra- 
tios I Trager et al. 2000b) and possibly even age (e.g. lThomas et al.1 
l200alNelan et al.ll2005i) . 

Unfortunately we do not currently have detailed predictions of 
velocity dispersions available for our model galaxies. However, it 
is well known that in the local universe, there is a tight correlation 
between stellar mass and velocity dispersion for early type galaxies. 
We therefore simply use an empirical conversion. 



logcr = 2.2 - 0.1023 



21.49 



'logM, - 1.584 
. 0.433 



(4) 



This relation is th e result of fitting to th e stellar mass and velocity 
dispersion data of Haring & Rid l l2004h and is consistent with the 
results of Lintott et al. (2006) obtained from SDSS. This parametri- 
sation of cr is a better approximation to the stellar velocity disper- 
sion than estimates that simply assume that it is proportional to the 
velocity dis persion of the dark m atter (sub)halo within the virial 
radius (e.g. ISpringel et aObOOlh . Early-type galaxies are baryon 
dominated so the stellar velocity dispersion can scale quite differ- 
ently from the (sub)halo velocity dispersion. 



2.8 Creating synthetic spectra 

We convolve the resulting star-formation histories (SFHs) of 
the m odel galaxies with t he m ulti-metallicity stellar SED mod- 
els of 'Bruzual & Chariot' ("2003, hereafter BC03), base d on the 
Padov al994 tBertelli et al.. ,1994) isochrones with a IChabrieil 
( l200lh initial mass function (IMF), to produce synthetic spectra, 
broad-band magnitudes and colours, absorption-line strengths, and 
mass-to-light ratios. 

The predicted magnitudes, line-strength indices, mass-to-light 
ratios, light-weighted stellar population parameters, and inferred 
SSP-equivalent stellar population parameters of the galaxy forma- 
tion models are computed in a 'post-processing' step by combining 
the SFHs with the stellar population models. In detail, the SFHs 
provide a mass in stars formed at each point in a grid of ages (t) and 
metallicities ([Z /H]) for each galaxy in the model at the redshift of 
interest (in this paper, 2 = 0.023; see below). The magnitudes in 
photometric band k are computed as 



Mk 



-2.5 log 



(5) 



where i, j are the age and metallicity bins, A4i.j is the mass formed 
at age i and metallicity j, I^ j is the intensity per unit mass in band 
k for the stellar population model with age i and metallicity j, and 
ZPfc is the magnitude zero-point for band k. We compute both dust- 
extinguished and unextinguished magnitudes. We model the impact 
of dust on the magnitudes using an approach similar to that pro- 
posed by Chariot & Fall ( 2000) as modified bv lDe Lucia & BlaizotI 
( I2OO7I) . We include two dust components, one due to 'cirrus' in the 
disc and one associated with the dense 'birth clouds' surrounding 
young star forming regions. The face-on optical depth of the 'cir- 
rus' dust in the V-band is given by 



Tv,0 C)C Tdust.O ^cold"lcold/(»'gas) , 



(6) 
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where Tdust.o is a free parameter, Z^oid is the metallicity of the 
cold gas, mcoid is the mass of cold gas in the disc, and rgas is the 
radius of cold gas in the disc, which we assume to be 1.5 times the 
scale radius of the stars. We then use a 'slab' model to compute th e 
inclination-dependent extinction (see lSomerville & Primackll999l) . 
Young stars (t < 10^ yr) are additionally enshrouded in a screen 
of dust with optical depth t-bc,v ~ ^'B C'v.o- We take^ec = 3. 
We assume a Galactic extinction curve ('Car demetal]|l989h for 
the cirrus component and a power-law extinctio n curve Ax oc 
(A/5 500 A)", with n = 0.7, for the birth clouds jCharlot & Fail 
I2OOOI) . 

Line-strength indices are computed similarly to the magni- 
tudes as 



EW = A 1 



L 



c 



(7) 



where A is the width of the index in A, F^'-' is the flux in the line 
centre for the stellar population model with age i and metallicity 
j and F^'"* is the flux in the continuum for that model, for indices 
measured in A, and as 



Mag = -2.5 log 



L 



(8) 



for indices measured in magnitudes ( IWorthevlll994h . We make no 
attempt to account for dust in the line-strength indices ( iMacArthuJ 
l2005t) . 



TFD08. If uncertainties in the stellar population parameters are de- 
sired, they are computed by taking the dispersion of stellar popula- 
tion parameters from 10000 Monte Carlo trials using the provided 
index errors as l-cr Gaussian deviates. We compute the stellar pop- 
ulation quantities for both the entire galaxy ('total') and its bulge 
component ('bulge'). Finally, we compute light-weighted stellar 
population parameters for each photometric band by combining the 
star formation histories with the SEDs. 



2.10 Coma Cluster model 

For the analysis in the current paper we are interested in compar- 
ing our models with real observations of early-type galaxies. The 
Coma cluster is a convenient target, as a large body of reasonable- 
to-excellent quality line-strength data is available (e.g., see the 
compilation in TFD08). We therefore generated twenty realiza- 
tions of a Coma cluster-sized halo, i.e., a halo with a circular 
velocity Vc ~ 1066 km s~^, corresponding to a virial mass of 
Mvir ~ 10^*'® M0, similar to if slightl y lower than recent es- 
timates of the m ass of the Coma cluster fcokas & Mamonll2003l : 
iKubo et al.ll2007l . and references therein). The smallest halos that 
we resolve in our merger trees have a virial mass of 10^" M0, so 
we select galaxies with M, > 10^ Mq, corresponding to a velocity 
dispersion cr = 57 km using Eq.|4] We select only 'early-type' 
galaxies (E and SO) according to their B-band bulge-to-total light 
ratios, as described in Section |23 
simulations to be z = 0.023 (e.g. 
redshift of the Coma cluster. 



We set the output r edshift of the 



Hudson et alj200ll) to match the 



2.9 Determining SSP-equivalent parameters 

As discussed in the Introduction, the most common approach for 
interpreting observed spectral line-strengths in terms of physical 
parameters that probe galaxy formation is to assume that the stellar 
population in a galaxy is a 'simple stellar population', i.e., that it 
is comprised of stars of a single age and metallicity. Under this as- 
sumption, one can use stellar population models to derive the SSP- 
equivalent age and metallicity that fits the observed line-strengths. 

To break the age-metallicity degeneracy, we must use both 
a metal line and a hydrogen line. For the analysis herein, we 
use the metal-line indices Mgb, Fe5270, and Fe5335, which are 
more sensitive to metallicity than to age, and a B aimer-line in- 
dex, which is more se nsitive to age than metallicity ( Worthev 1994; 
IWorthev & Ottavianiiri997h . Throughout the present paper we use 
the H/3 index, as this is nearly always combined with the previously 
mention ed metal-line indices to determine SSP-eq uivalent parame- 
ters (e.g.l Trager et al.ll2000al : iThomas et 'ZII2OO5'). Using a variety 
of other age and metallicity indicators is beyond the scope of this 
paper. 

We determine stellar population parameters directly using 
a non-linear least-squares fitting c ode based on the L evenberg- 
Marquardt algorithm described in iTrager et alj j2008l hereafter 
TFD08), in which stellar population models are linearly interpo- 
lated in (logt, [Z/H], [E/Fe]]f| on the fly to produce line strength 
indices (and magnitudes, colours and other parameters, but these 
are not typically used in the fitting process). For this purpose we 
use the BC03 stellar population models modified as described in 



^ Recall that the galaxy formation models have solar abundance ratios and 
therefore [E/Fe] = 0, although this is not true for the galaxy observations 
described later. 



3 SINGLE-STELLAR-POPULATION PROPERTIES OF 
GALAXIES IN HIERARICHAL MODELS 

In this section we seek to understand the properties of composite 
stellar populations when interpreted using single-stellar-population 
(SSP) equivalent par ameters. Thi s can be thought of as extending 
the work of Serra & TrageJ ( I2007D to more complicated - and hope- 
fully more realistic - star-formation histories. We first examine the 
SFHs and metallicity distributions of galaxies drawn from semi- 
analytic models, then examine their absorption-line strengths and 
inferred SSP-equivalent ages and compare these with the intrinsic 
physical timescales predicted by the semi-analytic models, such as 
luminosity- and mass-weighted ages, merger times, and birthrates. 

3.1 Star-formation liistories and metallicity distributions 

We begin by displaying the average SFHs and metallicity distri- 
butions of early-type galaxies in Coma cluster-sized halos in Fig- 
ure [T| In the upper panel we find clear evidence for 'downsizing' 
in the peak star-formation epoch of early-type galaxies drawn from 
the semi-analytic models: the more massive the galaxy, the earlier 
(on avera ge) the peak in the SFH (compare with the top panel of 
Fig. 1 in iDe Lu cia et al. 2006, who found similar results using a 
similar hierarchical galaxy formation model including AGN feed- 
back, although many of the details are different). The inset in the 
upper panel shows the SFHs of the galaxies over the last 2 Gyr, 
the time-scale to wh ich H/3 and the other Balmer lines are maxi- 
mally sensitive (e.g.. lSerra & Trager 2007, and Fig.|7]below); note 
that the lower the mass, the more likely the galaxy is to still have 
some ongoing star formation at these late times. In the current 
models however, downsizing is not very apparent at masses below 
K'h < 10^^ Mq, and galaxies with masses A/» ~ 10^^ M© have 
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Figure 1. Average star-formation histories (top) and metallicity distribu- 
tions (bottom) for early-type galaxies in 20 realizations of a Coma cluster- 
sized halo, binned by stellar mass. The inset in the top panel shows the re- 
cent low-level star formation in lower-mass model galaxies. Time in these 
panels is age or equivalently look-back time. 'Downsizing' is clear in these 
simulated galaxies: the age of the peak of the star-formation history de- 
creases with decreasing mass, although this decrease flattens out at the 
smallest masses. A mass-metallicity relation is also apparent in the average 
metallicity distributions, with a 'kink' and flattening in the peak metallici- 
ties around A/« ~ 10^^ Mq. 



the highest fraction of recent {t < 2 Gyr) star formation. In the 
lower panel we see clear evidence for a mass-metallicity relation 
in the early-type galaxies: low-mass galaxi es have lower metallic- 
ities than high-mass galaxies, as expected jFabeJll973h . Further- 
more it is clear that there is a 'kink' in the relation at a mass of 
M« ~ 10^^ M© and a flattening above this mass (note the spac- 
ing of the peak of th e metallicity distribu tions; compare also Fig. [8] 
below with Fig. 8 in lGallazzi et alj2005l) . 




log M./Mg 



Figure 4. Mass-age relations for model early-type galaxies drawn from 20 
realizations of a Coma cluster-sized halo. Top: mass-weighted stellar age as 
a function of stellar mass. The solid line is a Scr-clipped least-squares fit. 
The relation is quite shallow and possibly even flat below M, = IQI^Mq. 
Bottom: SSP-equivalent stellar population age as a function of stellar mass. 
The dashed line is a Scr-clipped fit to these points, and the solid line is the 
fit from panel (a). Note the large scatter in tssp> the lower appai'ent ages, 
and the strong apparent 'archaeological downsizing', which exaggerates the 
trend actually present in the mass-weighted ages. 



decrease with decreasing velocity dispersion, a phenomenon often 
referred to as 'arch aeological downsizing' (cf. iNelan et al ] |2005l: 
IXhomas et al.l2005l among others). The scatter also grows with de- 
creasing velocity dispersion, suggesting that the formation of low- 
mass galaxies is typically more extended than that of higher-mass 
galaxies (as we saw in Figure [T). The SSP-equivalent metallicities 
show a tight and strong mass-metallicity relation, as expected (see 
above), but the slope is steeper and the zeropoint is lower than the 
real Coma cluster data shown by the dashed lines; we discuss this 
point further below and in Sec. 14.41 



3.2 Line strengths and SSP analysis 

Using the procedure detailed in Section [2^ above, we have com- 
puted line-strength indices for simulated early-type galaxies in 
Coma cluster-sized halos. The results for H/3, Mgfe, and (Fe) are 
shown in Fig. |2] The left-hand panel shows that there is a weak 
trend in metallicity as a function of age, although with a large scat- 
ter. This is in the sense that lower-metallicity galaxies have younger 
ages, and originates mainly from the trends between galaxy mass 
and metallicity and galaxy mass and age in the models (low-mass 
galaxies are both younger and more metal-poor). The right-hand 
panel of Fig.[2]shows only that the enhancement ratios of the model 
galaxies are solar, i.e., that [E/Fe] = 0. This is not surprising, since 
we have used solar abundance SSP models to generate the synthetic 
spectra. 

Figure [3] shows the three derived stellar population parame- 
ters (age, metallicity, and [E/Fe]) plotted against the galaxy ve- 
locity dispersion. The SSP-equivalent ages of the model galaxies 



3.2.1 'Archaeological downsizing' and recent star formation 

Given the known sensitivity of B aimer lin es to recent s tar for- 
mation (recent star formation: RSF, e.g., iTrager et alj l2000bl; 
ISerra & Trageill2007t) . we ask whether the 'archaeological down- 
sizing' seen in these models is a reflection of true downsizing in 
the mass-weighted ages of the (model) galaxies or just results from 
small amounts of star formation occurring preferentially in smaller 
galaxies with extended star-formation histories. In Figure|4l we see 
that the latter case is likely to be the dominant signal in 'archaeolog- 
ical downsizing' in these models: low-mass galaxies have more re- 
cent star formation than high-mass galaxies, as expected from Fig- 
ure [T] resulting in younger SSP-equivalent ages. This can be seen 
even more clearly in the left-hand panel of Figure |5] which shows 
that the SSP-equivalent ages of galaxies are always younger than 
the mass-weighted ages, with the discrepancy (and scatter) grow- 
ing towards younger - and thus lower-mass, following Figure |4] 
- galaxies. We therefore suggest that 'archaeological downsizing' 
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Figure 2. Line-strength indices of simulated early-type galaxies from 20 realizations of a Coma cluster-sized halo. Left: The H/3-[MgFe] diagram is a 
convenient plot from which to read SSP-equivalent age (solid , roughly horizontal lines) and SSP-equi valent metallicity (dashed, roughly vertical lines), as 
these indices are at best only mildly affected by [E/Fe] 7^ jTrager et al.ll2000"3 ; lThomas et al ] |2004) . Right: The Mgb-(Fe) dia gram is a convenient plot 
from which to read [E/Fe]. The left-hand grid has [E/Fe] = (i.e., solar composition) and the right-hand grid has [E/Fe] = +0.3. Note that these model 
galaxies show no evidence for [E/Fc] 7^ 0, as expected from the chemical evolution prescription used in the simulation (see text). 
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Figure 3. Velocity dispersion-stellar population parameter relations for model early-type galaxies selected from an ensemble of 20 realizations of a Coma 
cluster-sized halo. The dashed lines are projections of the observed Coma cluster Z plane (for metall icities of [Z/H ] = —0.5, 0, +0.5 dex in the left panel and 
ages of 5, 10, 15 Gyr in the middle panel) and [E/Fe]-(T relation (right panel) on to the relations (see lTrager et al.l20 00b: TFD08). One can see 'archaeological 
downsizing' in the left panel, along with a steeper mass-metallicity relation than seen in the data (middle panel). Note that the point with high [E /Fe] in the 
right-hand panel is also the object with very young age and high velocity dispersion (a ^ 300km s~^) in the left-hand panel: it falls in a sligtly unusual 
region of line-strength space because of a burst of recent star formation and has a biased [E/Fe] as a result, even though intrinsically it has [E/Fe] = 0. 



will always appear stronger than any true (mass-weighted) down- 
sizing in the stellar populations due to the bias caused by recent 
star formation in low-mass galaxies. This is not to say that down- 
sizing does not happen in these models - it clearly does, as shown 
in the top panel of Figure |4]- rather th at 'archaeolog ical downsiz- 
ing' results (such as Nelan et al.ii20()5l ; lThomas et al.u2005) almost 
certainly overstate the true magnitude of downsizing, as suggested 
by TFD08. 

Further examining Figure |5] we find that the SSP-equivalent 
ages of model galaxies are in fact correlated with both the mass- 
and light-weighted ages of the galaxies, but the scatter is large in 
both cases (especially for the mass-weighted ages), as expected 
from the preceeding discussion. Simple linear fits to the points in 
Figure[5]give the following relations: 

logissp = 1.51 (±0.02) log(t)„,ass- 0.568 (±0.002) (9) 
logissp = 1.31 (±0.01) log(t)iight - 0.351 (±0.001) (10) 

for 1801 model galaxies, where the uncertainties are formal la er- 
rors assuming equal weighting of points in the least-squares fits and 
the RMS scatters are 0.08 and 0.05 dex, respectively. We note the 
following points. 



• As discussed above, RSF strongly biases the inferred ages, 
such that SSP-equivalent ages are younger than mass-weighted 
ages by more than 40 per cent on average and younger than light- 
weighted ages by roughly 25 per cent on average, thus greatly ex- 
aggerating the effect of downsizing on the stellar population ages. 

• The scatte r is asy mmetric in both diagrams. As discussed by 
ISerra & Tragerl j2007l) . galaxies that appear young in tssp are not 
guaranteed to be - and in these simulations, unlikely to be - young 
in a mass- or light-weighted sense. 

• The biases can be expected to be even worse for smaller sam- 
ples of galaxies. As an example, we consider the relations between 
SSP-equivalent and m ass- or light-weig hted age for a mock cata- 
logue representing the lMoore et al.l 1 120021) sample of Coma cluster 
galaxies (see Sec. |4.3| below). For this sample of 108 (model) galax- 
ies, the slopes of the relations are steeper [logfssp cx (1.89 ± 
0.25) log(f)niass andlogtssp oc (1.53 ±0.15) log{t)iight] and the 
scatters are larger (0.18 and 0.16 dex, respectively). 

If the SSP-equivalent age of a galaxy does not accurately rep- 
resent its mass- or light-weighted age, is there some other relevant 
time-scale of its star-formation history that it does reflect? Using 
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Figure 5. SSP-equivalent age as a function of mass-weighted (left-hand 
panel) and (V-band) light-weighted (right-hand panel) age for model early- 
type galaxies drawn from 20 realizations of a Coma cluster-sized halo. Solid 
lines are lines of equivalence for the plotted parameters and dashed lines 
are linear-least-sequares fits given in Eqs.|9](left-hand panel) and[Tol(right- 
hand panel). The SSP-equivalent and light-weighted ages are correlated, as 
are the SSP-equivalent and mass-weighted ages, but the scatter is very large, 
the slope is not unity, and the SSP-ages are biased low. The scatter is larger 
for the more physically interesting mass-weighted ages. These biases are 
due to residual star form ation in some fraction of the early-type galaxies 
(cf . ISerra & Tragej2007h . 
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Figure 6. Comparison of SSP-equivalent age with galaxy assembly time 
scales. Top panels: igsp as a function of the time of last merger (left-hand 
panel) and the time of last major merger (right-hand panel). Bottom panels: 
tgsp as a function of the time at which 50 per cent of stellar mass was 
formed (left-hand panel) and the time at which 80 per cent of stellar mass 
was formed (right-hand panel). 

our knowledge of the mass- and star-formation histories of the 
model galaxies, in Figure |6] we present tssp as a function of the 
times of the last merger (mass ratio greater than 1:10) and the last 
major merger (mass ratio greater than 1:5), and the times at which 
50 per cent and 80 per cent of the present-day stars were formed. 
The correlations of tssp with these quantities is even worse than 
with the mass- or light-weighted ages. In the case of the merger 
times, recent dry mergers leave tssp relatively unchanged, while 
minor but gas-rich mergers can seriously perturb tssp- The lack of 
a strong correlation of tssp with the times at which 50 and 80 per 
cent of the present-day stars formed is due to again to the presence 
of RSF. If the last small fraction of stellar mass is formed recently 
in a galaxy, this RSF results in a young tssp- 
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Figure 7. Correlation of SSP-equivalent age and 'birthrate' parameter, 
which is the star formation rate averaged over a relatively short timescale t, 
divided by the star formation rate over the whole lifetime of the galaxy. 
Top left-hand panel: biilhrate with r = 10 Myr (60.01 Gyr)- Bottom 
left-hand panel: birthrate with r = 100 Myr (60.1 Gyr)- Top right-hand 
panel: birthrate with r = 1.0 Gyr (bi.oGyr)- Bottom right-hand panel: 
birthrate with r = 1.0 Gyr (61,5 Gyr)- (Galaxies with no star forma- 
tion within the appropriate time-scale are not plotted.) The correlation of 
tggpwith bo.i Gyr. oGyr. and fei.sGyr highlights the sensitivity of the 
SSP-equivalent age to recent star formation. 

Rather than a specific event in the history of the galaxy, the 
SSP-equivalent age represents a rough measure of the birth-rate of 
stars over the last 1-1.5 Gyr. This can be seen in Figure|7] where we 
plot the average birth-rate bt - the fraction of stellar mass formed 
over some time with respect to the total stellar mass - for four 
time-scales: 0.01, 0.1, 1.0, and 1.5 Gyr. We see that the best cor- 
relation of tssp with bt is when 1 < t < 1.5 Gyr, suggesting 
that SSP-equivalent ages inferred from H/3 are maximally sensitive 
to RS F in the past Gyr (as shown from simpler, two-burst mod- 
els in ISerra & TrageJbool) . We therefore conclude that tssp is 
mostly measuring RSF over ~ 1-Gyr timescales, not mass- or light- 
weighted ages, not merger times and not the time at which some 
(large) fixed fraction of stellar mass has formed. 

3.2.2 Metallicities 

We turn briefly to the metallicities of the model galaxies. We 
present their mass-metallicity relations in Figure [8] both for the 
mass-weighted and SSP-equivalent metallicities. As expected from 
the bottom panel of Figure[T] a distinct flattening of the relations oc- 
curs at masses above 1ok (M«/M(;)) > 11-3. Th is is in rough accor- 
dance with the results of iGallazzi et al.l ( l2005h . who found a simi- 
lar effect in the r-band-luminosity-weighted metallicities of SDSS 
galaxies (see also Fig. 6 of SOS). Our (l/-band) light-weighted 
metallicities are slightly too low compared with the data; the differ- 
ence in the scatter is due to observational uncertainties in the SDSS 
data that we have not attempted to model. However, our intent here 
is not to make a detailed comparison with the SDSS data; rather, 
we will generate mock catalogues to compare with the Coma clus- 
ter galaxies in Section |44l below. 

We compare mass- and light-weighted metallicities with 
SSP-equivalent metallicities in Figure |9] Unlike in the case of 
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Figure 8. Mass-metallicity relations for model early-type galaxies drawn 
from 20 realizations of a Coma cluster-sized halo. Top: mass-weighted 
metallicity as a function of stellar mass. The solid line is a 3cr-clipped 
least-squares fit to the relation. There is clearly a very strong, tight mass- 
metallicity relation which flattens out at masses above M, 3x10" Mq. 
The long-dashed (red) line is th e mass— r-band - light-w eighted metallicity 
relation from the SDSS sample o f iGallazzi et allfeOOSl) and the dotted (red) 
lines are the Icr marginalised confidence bounds on this relation. Middle: 
V-band light-weighted metallicity as a function of stellar mass. Solid, long- 
dashed and dotted lines are the same as in the top panel. Note that the model 
light-weighted metallicities are slightly too low compared with the SDSS 
galaxies. Bottom: SSP-equivalent stellar population metallicity as a func- 
tion of stellar mass. Solid, long-dashed and dotted lines are again the same 
as in the top panel. The scatter is slightly larger, the slope is slightly flat- 
ter and the zero point is slightly lower than for the mass-mass-weighted 
metallicity relation. 
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Figure 10. A nearly face-on (left-hand panel) an d the edge-on (right-hand 
panel) view of the Z plane jTrager et alj|2000bh for model galaxies from 
the ensemble of 20 realizations of a Coma cluster-sized halo. The dashed 
lines in the left panel are projections of the observed Coma cluster Z 
plane (TFD08) onto the age-metalHcity plane for cr = 50, 150, 250, 
350 km s^^. The solid fine in the right panel is the edge-on view of the 
Z plane. At a given velocity dispersion, model metalHcities are typically 
too low with too steep of a slope, as seen in the middle panel of Fig. [5] 

the right-hand panel shows that the slope of the edge-on view of 
the Z plane of the model galaxies is too shallow compared with 
the observations. The former problem is due to the chemical evo- 
lution model used and assumptions about the IMF. Arrigoni et al. 
(in prep.) have shown that a more detailed, multi-element chemical 
evolution model that includes enrichment by Type la supemovae 
and a slight flattening of the IMF (to produce more high-mass stars 
and therefore more metals) can remove this discrepancy. The latter 
problem is either due to the wrong mass-metallicity slope (which 
seems unlikely from Fig. [8]l or to a different correlation between 
tssp and [Z/H]ssp in the data and models. We consider these pos- 
sibilities in more detail below. 




Figure 9. SSP-equivalent metallicity as a function of mass-weighted (left- 
hand panel) and (V-band) light-weighted (right-hand panel) metallicity for 
model early-type galaxies drawn from 20 realizations of a Coma cluster- 
sized halo. Solid lines are lines of equivalence for the plotted parame- 
ters. The SSP-equivalent metallicities correlate well with mass- and light- 
weighted metallicities, albeit with a small offset. 

tssp, [Z/H]ssp correlates almost perfectly with mass- and light- 
weighted metallicity (as shown in Serra & Trager 2007), with very 
small scatter in the latter case. The SSP-equivalent metallicity is 
therefore a very good tracer of the light- or even mass-metallicity 
of a galaxy, modulo a small offset towards lower metall i city. 

Finally, we turn to the Z plane of iTrager et al.l jlOOObh , a 
relation between velocity dispersion and SSP-equivalent age and 
metallicity. In Figure[TO]we present these quantities projected onto 
the Z plane of Coma cluster early-type galaxies determined in 
TFD08. The left-hand panel shows that the metallicities of the 
model galaxies are too low at a given velocity dispersion, while 



4 COMA CLUSTER: OBSERVATIONS AND MOCK 
CATALOGUES 

We now to turn to our second and third goals, which are interpreting 
the properties of local early-type cluster galaxies in the context of 
hierarchical formation models and testing the models themselves 
using the data. 

4.1 Colour-magnitude diagram 

To begin our confrontation between the models and data, we ask 
whether the current galaxy formation models can reproduce the 
colour-magnitude diagram (CMD) of galaxies - regardless of mor- 
phology or colour - in the Coma cluster. In Figure [TT] we show 
the magnitudes and colours of a galaxies drawn from a single 
Coma-sized halo overlaid on the CMD of Coma cluster galaxies 
drawn from the si xth data release of the Sloan D igital Sky Sur- 
vey (DR6 of SDSS'Adelman- McCarthv et al .'2008'). We divide the 
CMD into red a nd blue regions using the magnitude-dependent cut 
of iBaldrv et alj ( l2004f ). which provides a fairly clean separation of 
both the model and Coma cluster galaxies. We see that the slope 
of the red-sequence of model galaxies matches the Coma cluster 
galaxies very well. The model galaxies appear to be slightly too 
blue, which is likely due to the too-low metallicities of galaxies in 
the model. There are also not enough truly blue galaxies down to 
the limits of the model and the spectral data. However, this may 
be due to the fact that the model galaxies are selected down to a 
fixed stellar mass, not a given luminosity. As blue galaxies have 
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Figure 11. Observed and predicted colour-magnitude diagram (CMD) of 
the Coma cluster. Large black dots show galaxies with re dshifts from the 
Sloan Digital Sky Survey (SDSS) DR6 iAdelman-McCarth v et al. 20081). 
placing them at the distance of the Coma cluster; small grey points are other 
galaxies in the same field without DR6 redshifts, assuming that they are also 
in the cluster; and (pink) diamonds are model galaxies drawn from a sin- 
gle reaUzation of a Coma cluster-sized halo, with no regard to morphology 
or colo ur. The dot-dashed green line demarcates the blue-red galaxy divi- 
sion of lBaldrvetal]|2004 ; note the excellent agreement of the model and 
observed red sequence slope and good agreement of the zero-points. The 
dashed vertical line represents the magnitude limit of SDSS spectroscopy, 
which is close to the model-galaxy mass selection limit of 10^ Mq. 



low mass-to-light ratios, star-forming galaxies with masses below 
our mass hmit can have luminosities higher than red galaxies above 
our mass limit. The good (if not perfect) match between the model 
and observed CMDs give us confidence to proceed with attempting 
to match detailed stellar population properties. We note that the lu- 
minosity function predicted by the model is also in good agreement 
with the observed Coma luminosity function. 



4.2 Observational material 

To compare the hierarchical galaxy formation models with ob- 
served stellar population parameters of early-type Coma cluster 
galaxies, we use the data presented in TFD08, concentrating on 
three samples: the LRIS data of TFD08 (hereafter 'the LRIS sam- 
ple'), which has excellent quality (in terms of signal-to-noise) but 
only twelv e galaxies; the large, morphologically-selected, complete 
sample of iMoore et al.l l2002l hereafter 'the Moore sample'), of 
somewhat poorer quality than the LRIS sample but containing ten 
times as many galaxies; and the large, red-sequence-selected sam- 
ple of iNelan et al.l | |2005L hereafter 'the Nelan sample'), again of 
poorer quality than the LRIS sample but again much larger. Char- 
acteristics of these three samples are presented in Table [T] 

The line strengths of the LRIS sample galaxies have been 
adjusted to represent those taken through an aperture with radius 
equal to on e-fourth of the effe ctive radius (r — re/4) of each 
galaxy (see iTrager etal.ll2008l for more details on the method, 
and note that the apertures in that paper are denoted in diameters, 
not radii). This has been done in an attempt to fairly compare the 
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Figure 12. Top panel: Velocity-dispersion function (VDF) of model galax- 
ies for a suite of twenty Coma cluster- sized haloes. Bottom panel: VDF of 
observed Coma cluster galaxies in the lMoore etal] i2002l) samp le (black) 
and of the mock catalogue built to match the lMoore eTal j2002l) galaxies 
(red hashed). Note that the VDF of the full suite is missing galaxies at the 
extrema of the observed VDF, and so these are also missing from the mock 
catalogue (see text). 



'global' stellar populations of the Coma cluster galaxies with the 
simulated galaxies, which are unresolved and therefore should be 
considered to represent 'global' properties. The re/4 aperture does 
not perfectly represent the 'global' properties, but we show in the 
Appendix that stellar populations observed through apertures with 
radii r = re are typically negligibly older and 0.10 dex more metal- 
poor. The line strengths of the Moore and Nelan samples were 
taken through fiber diameters of 2'.'7 and 2'.'0; the typical 'global' 
ages and metallicities (i.e., taken through an aperture of radius re) 
for these apertures are again negligibly older and 0.10 dex more 
metal-poor. We will refer to these shifts from observed to 'global' 
parameters below. 



4.3 Mock catalogues 

To make the desired comparison in a sensible way, we construct 
mock catalogues of early-type galaxies selected from the simula- 
tions. These mock catalogues consist of model galaxies selected 
to be early-type galaxies by the criteria discussed in ii2.3l above or 
red-sequence galaxies (see below) and to have the same velocity- 
dispersion function (VDF) as the observed galaxies. A VDF is 
merely the histogram of (logarithmic) velocity dispersions log a 
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Table 1. Coma cluster observations 



Sample 








rr^OO (A) 
'^Fc5270 ^'^> 


n-200 (A) 


TFD08 


12 


0.032 


0.033 


0.037 


0.041 


Moore et al. (2002) 


121 


0.096 


0.095 


0.102 


0.163 


Nelan et al. (2005) 


100 


0.074 


0.080 


0.091 


0.105 



Cols. 3-6: Mean index uncertainties for a galaxy with a = 200 km s ^. 



with a specified bin width; in this paper we use a bin size of 
A log a — 0.05 dex. As much as is possible, the mock catalogue 
is constructed such that each bin in the VDF is sampled from the 
simulated galaxies with the appropriate velocity dispersions with- 
out replacement. Only when a VDF bin is less populated in the 
simulation than in the observations are galaxies selected from the 
simulation more than once (although these galaxies will have dif- 
ferent index strengths due to the assumed index errors; see below). 
If there are no simulated galaxies in a VDF bin populated by the 
observed galaxies, as can be the case at the highest (typically cD 
galaxies) and lowest velocity dispersions, no galaxies are added to 
the mock catalogue. As noted in Section l277] the mass limit of our 
simulations is M, > 10^ Mq, corresponding to o" > 57 km s~^. 
Galaxies with velocity dispersions lower than this limit are not in- 
cluded in the mock catalogues. An example of the result of our 
sel ection algorithm is g iven in Figure[T2]for the mock catalogue of 
the lMoore et all ( |2002[) sample. 

The Nelan sample is a colour- and magnitude-limited sam- 
ple designed to select red-sequence galaxies (Smith etal. 20(3)- 
Galaxies in this sample were selected to have Rc < 17 mag and 
colours redder than A{B ~ Rc) > -0.2, where A{B - Rc) is the 
offset in colour at fixed magnitude from the colour-magnitude rela- 
tion of the red sequence. The red-sequence colour-magnitude zero- 
point was a median fit to all galaxies in a given cluster with Rc < 
16 mag, assuming a red-sequence slope of — 0.06 mag mag^^. 
These selection criteria were also applied to the model galaxies, 
ignoring the brighter limit of the zero-point fitting, which in prac- 
tice makes little difference. One complication is that, for simplic- 
ity, B and Rc magnitudes were not computed directly in the mod- 
els, but were instead determine d from SDSS qri magnit udes us- 
ing the transformations given in lBlanton & RoweisI l l2007r) . Again, 
this makes little difference in practice. Note that these red-sequence 
galaxies can have any morphology and are not selected to be early- 
type galaxies. 

For each galaxy selected from the simulation to populate the 
VDF, the average errors in the line-strength indices of the observed 
galaxies in that bin are applied to the indices of that galaxy to de- 
termine its 'observed' indices and are labelled as its line-strength 
index errors. SSP-equivalent stellar population parameters are then 
determined from the 'observed' indices as described in Sec. 12. 81 



4.4 Mock catalogues: comparison with observations 

We begin by comparing the predictions of the models with the 
TFD08 LRIS sample. The predicted and observed line strengths are 
shown in Figure [T3] We see that the H/3 strengths of the observed 
and model galaxies are quite similar, but the [MgFe] strengths are 
too low. Figure [14] shows that this discrepancy is driven by the 
metallicity offset mentioned above - the metallicities of the model 
galaxies are too low by ~ 0.4 dex and have a too-steep velocity- 
dispersion-metallicity relation (right-hand panels). The difference 
in metallicities is not due entirely to aperture effects, as we have 
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Figure 13. Observed (black) and mock (red) line strengths of the TFD08 
sample. Model grids are the same as in Fig. [2] Note that the mock galaxies 
have a similar mean age (H/3) as the bulk of observed sample - although 
the observed scatter is larger - but the metallicities ( [MgFe] ) are shifted too 
low. 

shown in the Appendix that such effects can only change the metal- 
licity by ~ 0.1 dex. The left-hand panels of Figure [14] show how- 
ever that the observed and model (SSP-equivalent) ages are similar, 
although the models show a weak 'archaeological downsizing' ef- 
fect that is not clearly present in the observations. Note that TFD08 
claim that there is no 'archaeological downsizing' in this sample, 
but this is admittedly a very small sample. Examining the Z plane 
in Figure[T5] we find that, apart from the obvious metallicity offset, 
the current galaxy formation models produce a similar Z plane to 
that observed. However, the slight correlation between tssp and 
[Z/H]ssp in the models seen in the bottom left-hand panel of 
this figure may be the cause of the too-steep velocity-dispersion- 
metallicity relation seen in Figure [14] Again, the LRIS sample is 
very small, but the quality of the data suggests that a large survey 
of a cluster like Coma with very high-quality data would be very 
instructive to probe these models even further; we pursue this idea 
further below. 

We repeat this exercise with the Moore and Nelan samples in 
Figures 11611211 It is clear from these figures that the poorer data 
quality of the Moore and Nelan samples makes definitive conclu- 
sions difficult, but the same trends seen in the LRIS sample are 
generally present. Although the lowest-velocity dispersion model 
galaxies are missing due to the simulation mass limit (see above), 
there seems to be a lack of model galax ies with ages tssp < 4 Gyr 
compared with the lMoore et alj ( [2002h sample at a < 200 km 
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Figure 14. Velocity dispersion-stellar population relations for the observed 
(top panels) and mock (bottom panels) TFD08 sample (see Fig. [5). As in 
Fig. 1131 the age distribution of the mock galaxies is similar to the observed 
galaxies (left-hand panels), although a bit steeper as a function of velocity 
dispersion, while the metallicities of the mock galaxies are too low (right- 
hand panels). This difference in metallicity is not entirely due to aperture 
effects, which can contribute only ~ 0.1 dex of the difference (see text and 
the Appendix). The dashed lines are the projections of the observed Coma 
cluster Z plane (see Figs.lsland llOt . 






t (Gyr) [Z/H] 

Figure IS. The nearly face- on (left-hand panels ) and edge-on (right-hand 
panels) views of the Z plane iTrager et al.l2000bl) in the observed (top pan- 
els) and mock (bottom panels) TFD08 sample. Contours are 68 per cent 
confidence intervals; note that the age-metallicity coirelations are accentu- 
ated by but do not arise from the age-metallicity degeneracy. As in Figs. 1131 
and ll4l it is clear that the metallicity of the mock catalogue is too low, but 
the age distribution is similar. Again, the difference in metallicity is not en- 
tirely due to aperture effects (see text). The dashed lines are the projections 
of the observed Coma cluster Z plane (see Figs.lsland llOt . 
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Figure 16. As for Fig.[T3] but for the lMoore etal] \20Q% sample. 
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Figure 17. As for Fig.[T4] but for the lMoore etal] j2002l) sample. 



in the left-hand pa nels of Figure [TTl the effect is less severe for the 
iNelan et alj ( l2005t) sample (left-hand panels of Fig.|20|l. However, 
the errors in the observations are large enough to make a quantifi- 
cation of this resuh problematic. 

It is also clear from the right-hand panels of Figures [17] and 
[20] that the mass-metallicity relations of the mock catalogues are 
steeper than those of the real data, as the large number of galax- 
ies makes comparison of the relations easier than in Figure [14] al- 
though this result is also apparent there. We again see an offset in 
the metallicities in these figures larger than the likely aperture ef- 
fect. Using galaxy formation models with more detailed chemical 
evolution prescriptions by Arrigoni et al. (in prep.), we have found 
that the slope difference in the mass-metallicity relation is mostly 
corrected by the proper inclusion of SNe la in the chemical evo- 
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Figure 20. As for Fig.[T4] but for the lNelan et"ZI <2005l) sample. 

Figure 18. As for Fig. [15] but for the lMoore et"ai] j2002l) sample. 
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Figure 19. As for Fig.[l3] but for the lNelan et aP J2005h sample. 
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Figure 21. As for Fig.[T5] but for the lNelan et"ZI <2005l) sample. 



lution model, and that the metaUicity offset can be corrected by a 
shght flattening of the initial mass function. However, proper inclu- 
sion of this more detailed modelling here is beyond the scope of the 
current paper. The apparent smaller scatter in the metallicity ditri- 
butions of the simulated galaxies at low velocity dispersion in the 
lower right-hand panels of Figures [TT] and |20| is interesting, as the 
scatter in the data appears larger. It is likely that much of this scatter 
in the data is from observational uncertainties (e.g., poor emission- 
line correction and low signal-to-noise), but that may not explain 
all of the scatter, and the tightness of the model metallicities may 
be another problem in the current models. It would be i nteresting to 
attempt to reproduce the recent results o f ISmith et alj flOOSi) . who 
studied dwarf galaxies at high signal-to-noise in the Coma cluster, 
to pursue this further. 



We conclude from the comparison of these three surveys of 
Coma galaxies with their mock catalogues that the current semi- 
analytic models reproduce the present-day SSP-equivalent ages 
moderately well, although the agreement is not perfect. The model 
apparently does not have the correct star-formation histories for 
these objects, such that the 'archaeological downsizing' seen in 
model galaxies is stronger than in the best observations. However, 
the moderately good agreement between the SSP-equivalent ages 
of the mock catalogues and the data suggests that the star-formation 
quenching recipe implemented in the hierarchical galaxy formation 
model produces roughly the correct amount of recent star formation 
and is therefore a reasonable if imperfect representation of the pro- 
cess suppressing star formation in these galaxies. The model metal- 
licities are clearly problematic due in large part to known deficien- 
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Figure 22. The absorption-line strengths of two idealised nearby rich clus- 
ter galaxy surveys. Black points are model galaxies selected to have early- 
type morphology; red points are model galaxies selected to be on the red 
sequence (with the same hmits as the Nelan et al. 2005 sample). Model 
grids as in Fig. [2] Both mock catalogues have the same line-strength er- 
rors as the TFD08 sample; the morphology-selected sample has the same 
number of galaxies and VDF as the Moore et al. 12002, see Fig. [12] bottom 
panel) sample; and the red-sequence -colour-selected sa mple has the same 
number of galaxies and VDF as the iNelan et"al] j2005t) sample. The age, 
metallicity, and enhancement ratio distributions of these samples can be 
easily determined with such a large sample with such small uncertainties. 
Nearly all young, metal-rich galaxies in the red-sequence sample are dusty 
late-type galaxies. 



cies with the simple chemical evolution prescription used here. We 
expect to correct these deficiencies in upcoming models (Arrigoni 
et al., in prep.). 



4.5 Idealised nearby rich cluster surveys 

In the previous subsection we found that the currently available 
samples of Coma cluster galaxies are either too small or of too low 
quality to provide sharp tests of the hierarchical galaxy formation 
models. We now ask what could be learned from a spectroscopic 
surve y of a Coma-sized c lust er with roughly 10 galaxies (similar 
to the iMoore etalj|2002l and lNelan et alj|2005l samples) with the 
data quality of the LRIS (TFD08) sample. 

We create two different idealised large surveys: one selected 
solely by morphology to be early-type galaxies with any colour and 
one solely by colour to lie on the red sequence with any morphol- 
ogy (although both are mass-limited by the simulation itself). To be 
precise, we selected for the first sam ple early-ty p e galax ies from the 
simulations to have the VDF of the fMoore et alj j2002h sample and 
for the second th e red-sequence col our-magnitude selection func- 
tion and VDF of lSmith etal] ( |2004|) . In both cases we applied the 
line-strength uncertainties of the galaxy in the LRIS sample with 
velocity dispersion nearest to each model galaxy to the model line 
strengths. 

We plot the resulting line-strengths in Figure [22] and SSP- 
equivalent parameters as functions of velocity dispersion in Fig- 
ure [23] Note that the early-type galaxy sample has more young 
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Figure 23. The predicted velocity dispersion-stellar population relations 
for two idealised nearby rich cluster surveys. Top panel: morphology- 
selected sample. Points with large error bars in age (right panel) are those 
that sit very close to or off the bottom of the SSP model grids in Fig. 1221 
Bottom panel: red-sequence-colour-selected sample. As in Fig. 1221 the 
age, metallicity, and enhancement ratio distributions can be easily mea- 
sured from such surveys, but the 'archaeological downsizing' seen in mass- 
weighted age in Fig.|4]is difficult to see even in this sample. 
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Figure 24. As in Figs. ll8l and ll5l for two idealised nearby rich cluster sur- 
veys. Top panel: moiphology-selected sample. Bottom panel: red-sequence- 
colour-selected sample. In comparison with Fig. 1181 the error contours are 
small enough to see clear deviations from the plane, as well as any curvature 
that may be present. 
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(issp < 5) galaxies than the red-sequence sample; these are blue 
early-type galaxies, absent from the latter sample. Nearly all of the 
young, metal-rich galaxies in the red-sequence sample are dusty, 
late-type (spiral) galaxies in the simulation. On the other hand, it 
appears that there may be more low-mass, old but metal-poor galax- 
ies in the early-type galaxy sample than in the red-sequence sample. 
This is apparently a result of the magnitude limit of the latter sam- 
ple, which selects against the smallest and therefore metal-poorest 
galaxies; however, this is a small effect, as the magnitude limit of 
the red-sequence selection {Rc < 17) is close to the mass limit in 
our simulation for old galaxies at this redshift (Fig. lilt. 

It is clear from Figures I22U24I that a large, high-quality sur- 
vey of a Coma-sized cluster would provide a wealth of informa- 
tion on the precise distribution of stellar population parameters in 
early-type and/or red-sequence cluster galaxies, showing any devi- 
ations from or curvature of the Z plane and providing stringent tests 
of both hierarchical galaxy formation models and previous data. 
As a simple but interesting example, suc h a survey would answer 
the qu estion posed by TFD08 (see also ISanchez-Blazguez et al.l 
l2006al lbl): are the nearly constant ages of early-type Coma clus- 
ter galaxies just a selection effect resulting from observing only 12 
galaxies, or is this a general result for the entire cluster? 



5 CONCLUSIONS 

In this study we set out to use hierarchical galaxy formation mod- 
els to understand the properties of local galaxies and, conversely, to 
use those properties to test hierarchical galaxy formation models. 
To accomplish these goals, we make use of a state-of-the-art semi- 
analytic galaxy formation model that includes a well-calibrated 
parametrization of galaxy mergers, black hole growth and both 
'bright' and 'radio mode' feedback from these black holes. We 
coupled these cosmologically motivated models to stellar popu- 
lation models to produce synthetic spectra and computed line- 
strengths from these spectra. We then used the line strengths to 
determine 'single-stellar population-equivalent' (SSP-equivalent) 
ages and compositions in precisely the same way as for the ob- 
served spectra. 

We created 20 realisations of a simulated Coma cluster-sized 
halo for this study in order to compare our resul ts with recent stud- 
ies of earl y-type ga laxies in th e Coma cluster jMoore et al.ll200 j ; 
iNelan et a l. 2005; T rager et al ]|2008). We began by comparing the 
inferred SSP-equivalent ages and compositions of model early-type 
galaxies with the mass- and light-weighted ages and compositions 
of these model galaxies. We found that SSP-equivalent ages are (at 
best) biased tracers of mass- or light-weighted ages, and are always 
younger than than these 'true' ages. This is because hot, young stars 
contribute muc h more strongly per unit mass to the Balmer line s 
than old stars l lTrager et alj l20()0bl |2005| ; ISerra & Tragej 1200% . 
The bias of SSP-equivalent ages towards young ages implies that 
the 'archaeological' or 'fossil downsizing' seen in previous studies 
(e.g., Thomas et al. 2005; Nelan et al. 2005; Clemens et al. 2006) 
overstates the 'true' downsizing in the early-type galaxy popula- 
tion. Rather, we find that the SSP-equivalent age correlates most 
strongly with the fraction of stars formed within the past Gyr. 
On the other hand, the SSP-equivalent metallicity is an excellent 
tracer of the light-weighted metallicity and a very good tracer of 
the mass-weighted metallicity of a galaxy. This is because hot, 
young sta rs contribute little to the metal lines in a composite 
spectrum jTrager et al.ll2005l ; ISerra & Tragej[2007h . These conclu- 
sions strengthen previous suggestions about the behaviour of line- 



strength-derived ages and metallicities based on simple two-burst 
stellar population models but extend them to more complex and, 
presumably, more realistic star-formation histories from a fully- 
cosmological hierarchical galaxy formation model. 

We then extracted mock catalogues from the simulation to 
represent three recent li ne-strength surveys of early-type galax - 
ies in the Coma cluste r: iMoore et al] ( l2002h . iNelan et all ( l2005h . 
and T rager et alj j2008 ^. These mock catalogues were intended to 
help us both to interpret the observations and to test the galaxy 
formation model itself. Unfortunately, drawing significant conclu- 
sions about the star-formation histories of Coma cluster galaxies 
was not possible: the size of the high-quality sample ( TFD08) was 
too small, while the da ta quality of the large samples i jMoore et al.l 
I2OO2I ; iNelan et alj 12005,) was too poor We could however note 
some interesting problems with the model from discrepancies be- 
tween the mock catalogues and the observed samples. First, com- 
parison of the observed velocity-dispersion-metallicity relations 
with the models revealed that the models are currently producing 
neither the correct zero-point of these relations - the model galax- 
ies are too metal-poor - nor the correct slope. The former problem 
is likely to be due to the lack of SNe la nucleosynthesis in the cur- 
rent models, while the latter problem has two likely causes: (1) the 
star-formation histories of the model galaxies are not the same as 
real galaxies, such that the variation of fssp with velocity disper- 
sion is too steep in the models compared with the data (left-hand 
panels of Fig.ll4t. leading to an incorrect tssp-[Z/H]ssp relation 
(left-hand panels of Fig. |15t ; and (2) the assumed IMF may not 
be correct (Arrigoni et al., in prep.). That being said, the reasonable 
agreement between the SSP-equivalent ages of the mock catalogues 
and the observed samples suggests that, while imperfect, the star- 
formation quenching model implemented in the hierarchical galaxy 
formation models is a reasonable representation of the true process 
suppressing star formation in early-type galaxies. We suggest that 
a large, deep, high-quality survey of a cluster like Coma - easily 
achievable with the current generation of wide-field spectrographs 
on 8-10-metre class telescopes - would provide a sharp test for the 
models that is not yet possible with the current data sets. 

We have shown that the combination of hierarchical galaxy 
formation models with detailed stellar population models provides 
powerful tools both for understanding the stellar populations of 
early-type galaxies and for testing the hierarchical galaxy forma- 
tion models themselves. We look forward to better datasets and to 
more realistic models for both galaxy formation and stellar pop- 
ulations. As has been mentioned several times, we are currently 
working on one of these aspects, introducing a detailed treatment 
of multi-element chemical evolution into the semi-analytic galaxy 
formation models (Arrigoni et al., in prep.). These enhanced mod- 
els will allow us to probe the star-formation histories of galaxies 
- and test the galaxy formation models - much more robustly and 
thoroughly than previously possible. 
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Figure A l. Line-strength gradients in elliptical apertures derived from 
iGonz alez ( 1993), in the sense dlog //dlog(r/r'e); see text for more de- 
tails. Solid lines are error-weighted means of the index gradients. 
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APPENDIX A: APERTURE CORRECTIONS 

Unfortunately, the Coma cluster is not quite far enough away for 
its member galaxies to be completely contained in a reasonably- 
sized slit or fiber. A fair comparison of the observational data with 
the galaxy formation models requires corrections to be made for 
the fact that early-type galaxies have line-strength gradien ts (see, 
e.g., lDavies et al.ll993l : [Gonzalezll993l : lMehlert et alj200d , among 
many others). We must therefore estimate the effect of these gradi- 
ents on the stellar population parameters. 



Table Al. Mean line-strength gradients derived from lGonzalej jl993l) 



Index d ' 


og//(iIog(r/re) 


H/3 


0.019 ±0.051 


Mgb 


-0.044 ±0.018 


Fe5270 


-0.049 ±0.018 


Fe5335 


-0.045 ±0.022 


Table A2. Stellar population gradients 


Parameter 


rjA re 1''35 


logt 


0.01 ±0.05 0.01 ±0.05 


[Z/H] 


-0.10 ±0.03 -0.10 ±0.03 


[E/Fe] 


0.00 ±0.01 0.00 ±0.01 



To make the required corrections, we ap peal to th e high - 
quality line-strength gradient measurements of iGonzalej ( Il993h . 
Gonzalez measured gradients in 40 local early-type galaxies, typ- 
ically along both the major and minor axes. We have used his 
'elliptical-aperture-weighted' line-strengths (see his Chapter 5) and 
taken error-weighted averages, where possible, of the major and 
minor axis line strengths in apertures of radii r — re/16, re/8, 
re/4, and re/2 (if only one axis was available, we used that axis). 
We then fit linear relations of the form 

log J(r/re) = a + 61og(r/re), (Al) 

where /(r/re) is the line strength within an elliptical aperture 
with the fractional equivalent circular radius r/re. We define 
b = dlog/(r/re)/dlog(r/re) the line-strength gradient for each 
galaxy. To search for trends in these gradients as a function of, for 
example, velocity dispersion, we plot the results as a function of 
log(j(r — re/4) in Figure IaTI We find only minimal hints of 
trends with velocity dispersion and a very large scatter, so we use 
the error-weighted mean of the gradients from all galaxies for each 
index to estimate our aperture corrections. These means are tab- 
ulated in Table IaT] The final (logarithmic) line strengths at one 
effective radius - a reasonable guess at a 'global' value - are then 
taken to be 

log /(re) = log/(r) + (6) * log(r/re), (A2) 

where (b) is the mean gradient of index / and r is the radius of the 
original aperture: for the LRIS data, — re/4; for the Moore sample, 
r — l'.'35; and for the Nelan sample, r — I'.'O. The errors in the 
mean gradients are added in quadrature (after scaling by the loga- 
rithmic difference in the radii) to the original line strength errors. 

We estimate the mean shifts incurred in the stellar population 
parameters when scaling from the LRIS and Moore apertures to 
'global' r — re apertures. Unfortunately, we do not have measured 
effective radii for many of the galaxies in the Moore sample. We 
have however measured the line strengths and stellar population 
parameters of the LRIS galaxies in an equivalently-sized aperture 
as Moore et al. (200 2). The results a re presented in Table lA2l As 
expected (see, e.g.. lTrager et al.l l2000b), we find a slight difference 
in the metallicity, such that the galaxies are 0.10 ± 0.03 dex more 
metal-poor in the 'global' re aperture than in the smaller apertures. 
We find however negligible differences (less than 3 per cent in each 
parameter) between the ages and enhancement ratios in the larger 
and smaller apertures. 



